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Abstract 

We review results about the density of typical lattices in M n . They 
state that such density is of the order of 2~ n . We then obtain similar 
results for random packings in W 1 : after taking suitably a fraction v of a 
typical random packing a, the resulting packing r has density C (u) 2~™, 
with a reasonable C (y) . We obtain estimates on C (v) . 

1 Introduction 

The problem of filling the euclidean space K™ with equal non-overlapping spheres 
has a long and celebrated history. It was started by J. Kepler in 1610, who 
studied the hexagonal two-dimensional packing, and conjectured that the face- 
centered cubic three-dimensional one is the densest possible. Later on, I.Newton 
claimed that in R 3 there cannot be more than 12 non-intersecting balls of unit 
radius touching a given unit ball. In the case when the centers of spheres form a 
lattice (i.e., a discrete additive subgroup of K n ) the question of finding lattices 
with densities high enough was studied by C.F.Gauss. Starting from dimension 
9, the highest density of the lattice packing is still unknown. Naturally, even 
less is known when the ball arrangement is not of a lattice nature. The problem 
constitutes an essential part of 18-th Hilbert problem. Whether or not the 
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best density is achieved on the lattice arrangements is a major open question 
in the field. Quite recently T.C. Hales and S.P.Fcrguson seem to have proved 
that for n = 3 the largest density is that of the face-centered cubic lattice 



packing . The proof is very long and computer assisted. It is discussed in [Oe2|. 
Constructions of dense enough sphere packings are often based on rather subtle 
algebraic technique, sometimes using algebraic geometry and number theory 



methods (cf. |CSfl, p^B, pT[, pel]). 

In this paper we are not trying to solve any of the above problems. Instead, 
we address the question of the density of typical lattice packings and of typical 
random packings. Here, the word "typical" refers to the natural probability 
distributions, defined below. In essence, the probabilistic approach to the lat- 
tice packing problem is not new, and existing results on dense enough lattice 
packings are implicitly based on it (cf. Jm|, ||, |Sch| 1). So, the results about the 
typical lattices presented below, are rather of the review nature, though they 
are presented in a new way. The most important element in them is probably 
the concept of a typical lattice itself. The importance of these results for us is 
that they give the reference frame we need when we discuss the properties of 
random packings, the main topic of the present paper. 

As it is usual in the probability theory, the reasonable results can be obtained 
only in the limit when the number of degrees of freedom of the system goes 
to infinity. For the case of lattices that means passing to the limit over the 
dimension n — > oo. On the other hand, for random packings we have infinitely 
many degrees of freedom already in the (infinite volume) finite dimensional 
euclidean space R™, so no other limit is needed here. 

Below we fix some notation. In the next section we present our results 
concerning the properties of the typical lattice packings. The last section deals 
with the random packings. 

Let a be a locally finite subset of R™ . Let Vn C R™ denote a cube with the 
side 2N, centered at the origin. By erjv we denote the intersection a (~l Vn- Let 
d(CTTv) denote the minimal spacing between the points of gn : 



d(<rjv) = min \x — y\ , 

x,y£<y N :x^y 

while r{ajsi) = — d{a^) . We define the geometric (or sphere packing) density 
of the configuration a in the box Vn to be the number 

An (a) = wr ■ 

Here B (x, r) is a ball of radius r centered at the point x. 

The sphere packing density of an infinite volume configuration a is defined 

as 



A (a) — lim sup A^r (c) 

N^oo 
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This corresponds to the classical problem of packing equal non-overlapping 
spheres in R n . 

In the special case of lattice packings the set a is a discrete additive subgroup 
of R™ and the density can be rewritten as 

A( ff )- V(ff)n 



det (a) 



where v n (r) = r n v n (l) is the volume of the sphere of the radius r in R™, 
v n = v n (1) , and det (a) is the volume of the fundamental domain of a. 

The classical question then is about the largest possible value of A (a) for 
any locally finite subset a C R™ or for any lattice a C R™, and also how to 
construct the corresponding <r-s. 

In this paper we discuss a different problem. We are interested in the density 
of a typical (random) configuration cr, and that of a typical (random) lattice. 
Thus, our main results deal with the distribution of the random quantity A (<r) . 
Namely, we show that in the lattice case the geometric density of a typical lattice 
is of the order of 2~™ (see Theorem [l] below). For the random sphere packing 
the geometric density of a typical realization of a corresponding random field 
is zero. The problem becomes interesting if we allow ourselves to decimate the 
configuration and to throw away a fraction of it, containing 'bad' sites. Then it 
turns out that the remaining random configuration has the geometric density of 
the same order as in the lattice case (see Theorem |] below), for a proper choice 
of the bad set. 



2 Random lattices 

In this section we discuss the lattice case. In the first subsection we present the 
results concerning the case of very high dimension n, which turn into a simple 
relation in the limit n —* oo. In the following subsection we consider the case 
n = 2, when again the results can be expressed by simple relations, due to 
explicit computations. 

2.1 The case of large dimensions 

Let a is a lattice. Clearly, its geometric density remains the same after multi- 
plication by a scalar. Therefore, we can consider here only unimodular lattices, 
i.e., those with det (<j) = 1 . The space of unimodular lattices in R ra is natu- 
rally isomorphic to the symmetric space A n = SL n (R) / SL n (Z) . It is equipped 
with the Haar measure \i n and \i n (A„) is finite. We therefore can normalize 
it so that fi n (A n ) = 1. After such choice the density A (a) becomes a random 
variable. As the following theorem shows, in order to get a nontrivial limiting 
distribution for it, one should normalize it by the factor 2" : 

Theorem 1 Let 

F n (a?) = /Lt„ {a G A„ : 2"A (a) < x} 
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be the distribution function of the random variable 2™ A (cr) . Then 

Urn F n [x)=l- e-' x/2 . 

n — >oo 

This result follows from the following theorem of Schmidt: 

Theorem 2 Jffc^/ Let n > 13. Let S be a Borel set in W 1 such that Sn(-S) = 0. 
Suppose that vol (S) < n — 1 . TTien i/ie measure 

fi n {(7 6 A„ : (7 n 5 = 0} = e"™ 1 ^ (1 - i?„) , 

w/iere 

li^j < 6 f £) e 4vol ( 5 ) +vol(5) n - 1 n-"+ 1 e TOl ( 5 )+ n . 

■ 

Proof of Theorem |l|. To derive Theorem [j] from Theorem || take S — S (d) 
to be "one half of the ball" B (0, d), in such a way as to satisfy the conditions 
S n (-S) = and S U (-5) = B (0, d) \ 0. Then, for a lattice a the condition 
cr (~1 S — is equivalent to the condition that d (cr) > d, so 

H n {a £ A„ : d (a) > d } = fj, n {a G A„ : a n S (d) = 0} 
= e -™KS(d)) (1 _ Rn) . 

Let us put a; (cr) = v n d(a) n , then 2"A (cr) = a; (cr) , since det (cr) = 1. Put 
d = ( f- J . Then, by definition 

F n (x) = fl n {cr G A„ : 2™ A (cr) < x} 
= fJ>n W eA„:i (a) < x} 
= 1 - /!„ {cr G A„ : d (cr) > d} 
= 1 - (i n {cr G A n : cr n S (d) = 0} . 

Noting that vol (S (d)) = f , we have 

F„ (x) = 1 - /i„ {cr G A n : a n 5 (d) = 0} - 1 - e - vol( - s ^ (1 - R n , x ) 
= 1 - e"*/ 2 (1 - i?„, x ) , 

with 

|iin,x| < 6 (~) ^ e 4 ^ 2 + n-"+ 1 e^ 2 +" - 

as n — > oo, which proves our statement. 
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2.2 Two-dimensional case 

For lattices in R 2 the information is much more precise. 

Consider lattice packings in R 2 , i.e., lattices in C. Instead of fixing their 
determinant, let us use the freedom of scalar multiplication and rotation to fix 
the shortest basis vector to be 1 . Then we can take the other basis vector z in 
the modular domain 



, 1 1 



The corresponding lattice will be denoted by L z . Its packing radius is i. If 
z = x + iy, then the volume of the fundamental parallelepiped of L z is y. 
Therefore, its sphere packing density A (L z ) equals 

The space of two-dimensional lattices A 2 = SL 2 (R) / SL 2 (Z) is canonically 
homeomorphic to what we get from F by identifying every boundary point (x, y) 
with the boundary point (—a;, y) . 

We will use the well-known fact that the probability Haar measure on A 2 in 
the (x,y) coordinates is given by dfi = f-^pr 2 . The sphere packing density A 
then becomes a random variable. 

Theorem 3 The random variable A has a density pa (x), given by 



PA (x) 
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for < x < f 



4 ' 



i-V 1 "^) 2 ) for$<x<j-, 
forx>^ 2 . 



In particular, the mean value of A equals | log 3, its variance equals — ^ (l°g3) 2 ^ 

and the maximal possible sphere packing density equals -j— . 

Proof. The lattice L z has density A if and only if y — y (A) = j^. The 
parameter z = x + iy is therefore uniquely determined by the pair (x, A) . So it 
can be taken as the coordinate system on F. We have dy (A) = — f ^ ■ The 
measure /i can be rewritten in coordinates (x, A) as 

dfi= 3_dx^ = 3_dx^) = ^ xdA 

7T y 7T (^) 7r 

Let Fb — Ffl {y = b} . The density pa (a) is then equal to ^|mcs (F y ^) , where 
mes(-) stands for the one-dimensional Lebesgue measure. For y (a) > 1, i.e., 
for a < j, we have mes (i^( a )) = 1; f° r 2/ ( a ) < "^7 i- e -7 f° r a > ^757 we 
have mes (F y ^) = 0. For ^ < y(a) < 1 the total length mes (F y ^) equals 

2 — \Jl — y (a) 2 ^j = 1 — 2^/ 1 — (j^) 2 , which proves the first statement. The 
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statement about the maximal value of density follows immediately. The mean 
value is given by the elementary integral 



and a similar equality defines the variance. ■ 

3 Random packings 

In this section instead of lattices we will consider configurations a of the point 
random fields, and we will try to solve for them the same questions we were 
discussing above for lattices. 

By a point random field in R™ we mean a probability measure P on the set 
S of all countable locally finite subsets of R". The simplest example of such 
measure is the Poisson random field. To define it we first introduce for every 
V C R" the notation Sy for the set of all locally finite subsets of V } and for 
every a G <S we denote by ay G Sy the intersection a (~l V. We denote by \ay\ 
the cardinality of the set ay ■ 

A random field P A is called a Poisson random field with intensity A > if 
and only if 

i) for every two disjoint subsets V, W C R™, V D W = 0, the random configura- 
tions ay, aw are independent; 

ii) for every finite subset V C R™ the conditional distribution of ay under the 
condition that \ay\ = m is just the Lebesgue measure on V m , normalized by 
the factor , 1( L m , while the probability of the event | ay | = m is given by 



(We recall that the independence property means that for every two events A C 
Sy , B C S w we have P {a : ay G A, a w E B} = P {a : ay e A} P {a : a w G B} 
provided V, W are disjoint.) In what follows we will denote by P the Poisson 
random field with intensity A = 1, and we will omit A from our notation, as well 
as the adjective "Poisson" . 

More general random fields can be treated by the methods presented below. These 
are called Gibbs random fields corresponding to the interactions U. Here U — U (x, y) 
is some given function interpreted as the strength of interaction between two particles 
situated at locations 1/1/6 R n . A random field ¥ u is called Gibbs random field with 
interaction U if and only if its conditional distribution inside the finite box V C R™, 
given the configuration <7Rn\y outside it, has density with respect to the measure P x 
proportional to the "Gibbs factor" 




Avol(V) ( Avol (^)) 



m 



F{a : \<7y \ = m} = e 





> . 
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(For a general function U the Gibbs random field lacks the independence property 
i) above.) To make it into a well-defined object, one has to put some restrictions on 
the function U, The most general one is called superstability condition. We will not 
formulate it, but just note that it is satisfied if U is repulsive, which just means that 
U > 0. For further details see ]r| or 0. The Poisson fields above correspond to zero 
interaction. 

Let now a be our point random field. Then it is easy to see that for every 
e > 

P (A N (a) > e) -> as N -> oo. 

Since our goal is to construct subsets with positive (and even as big as possible) 
geometric density, this is unsatisfactory. To save the situation we allow to 
decimate the configuration a, by considering only "a fraction i/" of our collection 
of points <7, where v £ (0, 1) . To this end we introduce the quantity 

Aat,„ (cr) = max A N (t) , (1) 

tCctjv: 

H CTJ v H <7« | 1/2+e < M < H CTJV | + 1 <t n | 1/2+e 

where e is some fixed small number. This fraction is so designed that for any r 
satisfying the restrictions in (Q) we have rj^j-r — > v as N — * oo for a typical. As 
the following statement shows, this improves the situation. 

Theorem 4 With F -probability 1 the limit 

T> (y) — lim inf Ajy „ (cr) 



i) Let us introduce the function V\ (d) = 1 — \v n (d) . Then for every d > 



does not depend on a, for every is, < v < 1 

i) Let ms introduce th 
we have the lower bound 

V{v l {d))>2- n v n {d)y 1 {d), (2) 

provided of course that V\ (d) > 0. In particular, for d = d n , where d n satisfies 
v n {d n ) — 1, we <;ef v\ (d n ) — s, and so 

m) Introducing the function V2 (d) = 1 — (rf) e ^ 2Vn ^ d \ we have for every 
d>0 

V{v 2 {d))<2- n v n {d)v 2 (d). (3) 

ra) Introducing the function 1/3 (d) = 1 — (d) — (d) e~"™( d ), we /lave 
/or every d > 

X> (^ 3 (d)) > 2- n v n (d) ^3 (d) . (4) 
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Figure 1: Upper (3) and lower (2,4) estimates on V (v) 



Note 1. The relations (|])-(||) give implicit bounds on the function V [v) . 
The relation (J2J) can easily be rewritten in explicit form: T> [y) > 2~ n+1 v (1 — v) . 
The remaining relations can not be rewritten so easily. 

Note 2. The relation (|J) is an improvement of (||), see Fig. 1. 

Note 3. Because of the properties of the function V2 (d) the relation (||) 
gives an upper estimate on V (y) only for the values v > 1 — 4e _1 w 0.908 . 

Note 4. Further improvements of the relations (||)-((|) can also be obtained, 
see the proof of the theorem below. 

Note 5. The above theorem tells us that for almost every configuration a 
and for every N large enough we can find a subconfiguration a' N C <jn with 
about v |ctjv| points in it and with the geometric density close to T> (y) . However, 
the configurations <j' n might not converge as N —* oo. 

3.1 Proof of the Theorem 4 

We impose periodic boundary conditions in the box Vn, which wraps it into a 
torus. For a finite subset of points a C Vn, a point x £ Vn and a number d > 
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wc define 



m (x, d, a) — \B (x, d) n {a \ x}| . 

This is just the number of points in a, different from x, which are at a distance 
not bigger than d from x. Let us also introduce the measure 8 a on Vjv by putting 
a ^-measure at each point of a : 

The relevant quantity to look on is now the sum 



M (d, a) = — [ m (x, d, a) S a (dx) . 
2 Jvn 



M (d, a) is the number of pairs of points in a at a distance not bigger than d. 
The reason we introduce the quantity M (d, a) is the following. Let us consider 
the graph Gd (<r) with vertices at the points of a and with edges connecting any 
two vertices at distance < d, then M (d, a) is precisely the number of edges in 
Gd (<t) . Clearly, in the case M (d, a) = M, we can find a subset r C a with 

M = h-m, 

such that d (r) > d, so the graph Gd (r) has no edges and therefore 

Aiv(r)> ^ . (5) 

Irl 

Since evidently for v = - — r we have 

&n,u W) > Ajv (r) , (6) 

the rhs of the estimate (||) provides the lower estimate for the quantity we are 
interested in. 

We introduce now other important quantities, which allow us to obtain better 
estimates on the number of points in the subset r C a with the property that 
d (t) > d. Let Mi (d, a) be the number of isolated edges of the graph Gd (f) . 
Then for every r with d (r) > d we have 

|t| < |ff| - Aft. (7) 

Likewise, we introduce the quantity M3 (d, a) to be the number of subgraphs of 
Gd (c) which are maximal connected components of Gd (c) with 3 vertices and 
3 edges (i.e. just isolated triangles), and Mi (d, a) to be the number of maximal 
connected components of Gd (cr) with 3 vertices and 2 edges. Then we have for 
all r with d(r) > d that 

Irl < Id - Afi - Af 2 - 2M 3 . 
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On the other hand, because of the Theorem || below, we can claim the 
existence of the subset tCu with d (r) > d and such that 

\t\ > \a\ - Afi -\{M- Mi) (8) 



(see estimate (|24|) below), and even 

3 

1-7-1 > \a\ - Mi - M 2 - 2M 3 - - (M - Mi - 2M 2 - 3M 3 ) 

5 



(see estimate ( |25| ) below). One can proceed further with such estimates, in- 
troducing the quantities Mi with higher is; any information we have on the 
behavior of the random variables Mi (d, a) provides us with some answer to the 
question we are interested in. 

We start by the study of the random variable M (d, on) . 

Lemma 5 The mean value 

E(M(d,a N )) = ^-(2N) n . 
Proof. We first use the identity 

E(M(d,a)) = E [E (M (d, a) | {|er| = k})} . 

The conditional expectation E (M (d, a) \ {\a\ = k}) can be computed in the 
following way. Denote by \d (x, y) the function 



Xd (x,y) = 



1 if |x — y\ < d, 
if | x — y | > d, 



x,y G Vn- Let £i,...,£fc be k independent random points in Vn, distributed 
uniformly according to the Lebesgue measure on Vn ■ Then 



E{M{d,a)\{\a\ = k})=E\ £ X d&,Si) 

. l<i<j<k 



k(k-l) Vn (d) 

2 (27V)™ 



Hence 



E (M (d, a)) = £*(*-!) ^!!e-(-r , «*M W 

' fc>i 



In the same way one can compute the expectation E (Mi (d, &n)) • A particle 
x G (7/v contributes to Mi (d, cat) if there is another particle y G on with 
|rc — y| < d, and there are no other particles in the union B (x, d) U B (y, d) . 
Since vol (B (x, d) U B (y, d)) G [v n (d) , 2v n (d)] , we obtain 

V -^- (2N) n e~ 2v ^ < E (Mi (d,o N )) < (2N) n e~ v '^ . 
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Before dealing with the variance, we will make a slight generalization. Namely, 
let (f>(x,a) be a 'local observable', that is, a function which depends only on 

the intersection B (x, ctj n {a \ x} , for some d = d{4>) . We introduce a random 

variable $ (a) by 



$(<7)= / (f>(x,* N ) 5 a (dx). (9) 

JVn 

For example, the choice <j> (x, a) — m (x, d, a) corresponds to $ (a) — M (d, a) . 
In case 

!1, if m (x, d, a) = 1, and for the unique 
y e <7 with | a: — y| < d, m (y, d, a) = 1, 
in all other cases 

we obtain <!> (a) = Mi (d, a) . We will impose the following restriction on cf) : 
there exists a constant c > 0, such that for all x,y £ R n , all a C M" 

|^(i,<TU»)-^(a!,ff)|<c. (10) 

It clearly holds in the above examples, with c = 1. 

Lemma 6 T/ie variance D($(tTjv)) satisfies the estimate 

D($(fTjv)) < C(2N) n , 

with C = C((j>). 

Proof. We start with the identity 

D(0=B[E(C|r?)]+E[D(^|r ? )], (11) 

valid for any two random variables. We are going to apply it for the case of 
£ = $ (er) , while n will be the restriction of a to certain subsets K C Vn, which 
are called "corridors" and which are defined as follows. Consider the partition 
II of the box Vjv ; formed by the cubic subvolume Vi together with all its shifts 
by the vectors of the lattice 2lZ n . The number I is chosen to be equal to d + D, 



where the number D has to be of the order of d : D e 



dk(n) ,2dk(n) 



with 



the integer k (n) depending only on the dimension n. For such a partition to 
exist, we need the ratio to be an integer; clearly, such a choice of the 

d+D ° ° 

number D is possible. The union of the corridors K\ is now defined as the 
(^-neighborhood of the union of the boundaries of all the boxes of II. In other 
words, it consists of the strips of the width 2d, which are parallel to the various 
coordinate planes, with the spacing between the two consecutive one to be equal 
to 2D. (The subscript in the notation K\ is needed due to the fact that later we 
will have to consider other corridors.) The reason to introduce the corridors K\ 
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is that the random variable $ (a) , conditioned by the value ij of the restriction 
<j\k\i turns into the sum of independent random variables. 

Let us start with the second term of (pd|), E [D (£\r))] . We first estimate the 
contribution to the inner variance, coming from a single box Vi of II. That is, we 
need to consider the box VJ_j, the Poisson random field a inside it, together with 
the fixed configuration r\K(yi) m the corridor K (Vi) — Vi\V t _^. The random 
variable we should study, is the sum 

§ Vl (c I Vk(v,)) = 4>{x, a) S a (dx) + [4> (x, a U 57) - <f> (x, a)] 8 a (dx) . 

(12) 

To estimate its variance we will use the evident relation: 

D (C + x) = © (C) + © ( X ) + 2 (E (Cx) - E (C) E ( X )) (13) 



with £ to be the first term in (]12j), and x ~ the second. The first variance is just 
a constant: 



R = R (Z, (j>) = D ^ ^ (a:, a) ^ (dx) J . 

To estimate the covariance term, corresponding to E (£x) — E (£) E (x) in (|l3|), 
we rewrite the difference [<p (x, a U 77) — (a;, <r)] as 

n 

^ (as, o U ??) - (x, a) = ^2 i 1 ( x ' a u Vi-i,Vi) > 

i=l 

where 

■0 (a;, k, y) = (/) (x, k U y) - <p (x, k) . 

Here we use some enumeration ^(v,) = {j/ii 2/2, — , J/n} , while rji = {yi,y 2 , —,yi} ■ 
Then 

E(Cx)-E(C)E( X ) = 

r)S a (dx)l ip(x,aUr]i-i,yi)5 rr (dx) ] - (14) 



- E yJ^ ^(x,a) S a (dx)jE' 7 ^ ^(x,aUr) i - 1 ,y i )S <T (dx)J}. 

(The symbol E CT means taking expectation in <r; 77 is here a fixed parameter.) 
Note that 

ip (x,a U r)i-i,yi) = <p(x,a U 77,-1 U j/j) - <?!> (a;, cr U 77^—1) = 
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unless x £ B (vi , dj , in which case 

ip(x,aUrn-i,yi) = ip { x , ff B ( yi ,2d) u? 7>-i;^) 



and 



Therefore 






-E ff 




= E' 7 


(X, 


-E ff 





(x, c) <5 CT (dx) J ^ tp(x,aU ifr-i, y<) c5 CT (dx) I - 
(x, cr) <5 CT (cZx)^ E 17 (-0 (x, cr U T]i-i,yi) S a (dx)) 



Because of (pX)[), the random variable \tp (x, a U T/i—i, = 
= |0 (a;, c U 77i_i U — (^>(x, <T U 7jj_i)| is uniformly bounded by c (</>), while 

^(^B^d)) ^ C (^) °B( W ,2d) ' S0 E<T (^( aj '°'B(|f„2j)) ) iS bounded aS 

well. Therefore we can continue in (|lj) by 

E(Cx)-E(C)lE(x)<^i|»7K(v,)|. 



where R\ = R\ (<j>) , thus having 

O a {4> Vl (a | VK(Vl))) <Ro + Rl \VK( VI )\ + 



(15) 



it, yJ Vt ^ i x ' a B( yi ,2d) [J 'ni-i,yi) s a (dx)j\ 



Two random variables: A t = (j Vl ^ ( x > a s(y 2d) ^ Vi—li Vi) $a (dx)j and 
Aj = {jv l ( x > a B(y -2d) ^ Vj-uVj^j ^a(dx)J - are independent, provided 
\yi -yj\ > 2d. Hence, applying again @, 

pCT (E ^ ^ (^(^d) Um-uVi) S* (dx)^ = B a (l>0 ( 16 ) 

n 

= ^B CT A 4 + 2 £ [E(A i A j )-E(A i )E(A j )]. 



13 



The first term in @ can be estimated from above by 

R2(<P)\VK(v,)\, (17) 

while the second term - by 

Rati) f m(y,2d,v) S v (dy). (18) 
Jk(vo v ' 

Putting together ([TJ), ©, (0) and (Q, we get 

0((j> Vl {(T\m(v l ))) < Ro + R*\vk(k)\ + 



R 3 m 

IK(Vi) 



2d, 77) S n (dy). 



Hence the total variance 



yen 



< (t^—) R o + Ri [ S„ {dy) +R 3 [ m (y, 2d, a) 5 a (dy) 
\d + DJ Jki JKi k ' 

< ( . ) R0 + R4 [ S a (dx)+R 3 I m(x,2d,a) S„ (dx) , 
\d + Dj J Vn Jv N v ' 

and the expectation EJB(£|?7)] is indeed less than G\ (2N) n for some Ci, ac- 
cording to the Lemma g|. 

Going now to the first term, D [E (£\r))] , in ([Tl]), we will slightly abuse the no- 
tation, denoting by a the Poisson field outside the corridor K\, with r\ denoting 
the Poisson field on K 1 . We have 

E CT ($(aU77)|?7) = (^—] E ff ( / cf> (x, a) S a (dx) ) + 



E" Qf (a:, a U 77) - (a, a)} 6 aUv (dx)j 



The first expectation is just a constant: 

R b (0) = E uf (x, a) <5 CT (dx)^ , 

and so does not contribute to the variance. Therefore we need to study the 
variance of the random variable 



* 1 (r?) = E ff ( f [<j>(x,<rU v) - <f> (x, a)] 5 aUrj (dx) 
\Jv N 



(19) 



14 



We want to argue that this function can be represented similarly to (^): 

a 1 fa) = I ^(v,v) s n (dy), (20) 



K x 



with the function tp 1 (x,rf) satisfying the analog of relation jlCj). To show this 
let us first fix an enumeration rjK t — {yi, yi, Vn} for every configuration r]K 1 ■ 
(This even can be done in a measurable way, but since this enumeration will 
not be needed in the final analysis, we will not elaborate on this point.) Define 
Vi = {yi,V2, —,yi} and rewrite @ as 

$ x (77) = V E CT [ / 4> (x, a U Tji) 6^ {dx) - 
<p(x,a\Jr]i-x)5 a uni-i(dx)]. 

JV N 

We define now the function 

4> x (y,rf) = W[ct>(y,aUr h - 1 Uy) + (21) 
(^J 4> (x, a U ?7i_i U y) - t/> (x, a U S crUm _ 1 (dx) ] . 

It is almost what we need. Indeed, the relation (^) is straightforward, while 
the estimate ([To]) follows from the same estimate for </>. The only drawback is 
that the function 1 (y, rf) depends not only on rj, but also on the ordering we 
choose. However, if we now pass from (y, rf) to its symmetrization 1 (y,rj), 
obtained by averaging over all possible orderings on 77, we get what we want. 

The net result of the above discussion is that the estimate of the variance 
of the random variable $ (crjv) is reduced with the help of the identity (|ll| ) to 
the estimate of the variance of the random variable $ 1 (rj) . The key advantage 
of passing from $ to (& 1 is that the latter random variable depends only on 
restriction ij of a to "(n — l)-dimcnsional" subset K\ of K n . For example, in case 
n = 1 the set K\ splits into union of disjoint segments, and the random variable 
$ 1 (rf) is then the sum of independent random variables, corresponding to these 
segments. Hence, the variance of this sum is just the sum of the variances. For 
general values of n we have to repeat the above scheme for (n — 1) more times, 
reducing after each step the "dimension" of the corridors by one. 

So we proceed as follows. To estimate the variance of (rf) we again apply 
the identity jll|), with the following choices: £ = <f> 1 (cr|x 1 ), while 77 is the 
restriction of a to the "(n — 2)-dimensional" corridor K2 C K\, defined by the 
following construction. Let wi C Vi be the subset of 2n points of the cube 
Vi, consisting of the centers of all its (n — l)-dimensional faces. (The previous 
subset u>o C Vi contained just one point: the center of Vi itself.) Then 

K2 = K,\ |J (V^ + b). 

be{ivi+2lZ d } 
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(Again, for n = 2 the set K2 splits into disjoint union.) The key feature of thus 
defined corridor is that the random variable $ x (oi^i) conditioned by the value 
of the restriction <t\k 2 splits into the sum of independent random variables, so 
we can repeat the above arguments. The rest of the proof is just this repetition 
and is omitted. ■ 

Let £ be any random variable. The Chebychev inequality claims that for any 
e > 

P{|£-£(e)l>e}<^- 
Applying it to M (d, o\/v) we get the following 
Proposition 7 Let e > and S > 1/2 be fixed. Then 
\ M(d,a N )-^(2N) n 



((27V)T 



> e 



as N 



00. 



End of the proof of Theorem |l|. Let us fix some value of the distance 
d. Let N be large enough, e be fixed, and er/v be a typical configuration in Vn- 
Then it follows from the last proposition that with probability P going to 1 as 
N — > 00 the configuration ctjv has the properties: 



M(d,a N )-^(2N) n 



and 



((2A0 n ) 2/3 
Wn\ - (2N) r 



< e 



((2Nf) 



2/3 



< e 



(22) 



(23) 



(again by Chebychev). 

Plugging in these data into estimates (^), (@), we find that for v = 1 — ^v n (d) 



Ajv,!/ (c) > v n 



v n (d)j +o(iV- 1 / 4 ) 



uniformly over configurations ct/v, satisfying (|2^), (|23j), which proves (g). The 
proof of the rest of the Theorem follows by similar arguments, applied to rela- 
tions <§, (|). ■ 
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3.2 Vertex covering number of a graph 



In this subsection we prove a result from the graph theory, which was used 
above. The Theorem || below contains in fact a result slightly stronger than 
what was needed. 

Let G = (V (G) , E (G)) be a finite connected graph without loops and multi- 
ple edges. Let v (G) — \V (G) \ be the number of its vertices and e (G) = \E (G)\ 
the number of its edges. A set A C V (G) is called a covering vertex set if any 
vertex of G is a neighbor of a vertex of A. Another way to look at this property 
of A is to say that if we cross out all the vertices of A together with all edges 
having at least one end in A, no edges are left. The vertex covering number 
a (G) is defined as the smallest number of elements in a vertex covering set. 

Theorem 8 Let G be as above, with e (G) > 0, then 

a{&) < , 

and, in particular, for e (G) even 



Therefore, if e (G) > 2, then 



and if e (G) > 4, then 



Also, 



a{G)< l -e{G) 



a(G)<~e(G), (24) 



«(G)<jje(G). (25) 



a (G) f I + S f° rv odd ' 

max —77^ = S , i (26) 
G,v(G)=v e(G) y i + _L_^. for v even, 

where the maximum is taken over all graphs G of the described type such that 
v(G)=v. 

We say that a vertex is an end vertex if there is exactly one edge adjacent 
to it. 



Lemma 9 Let G be as above, with v (G) > 3. Then there exists a non-end 
vertex such that if we cross it out together with all its edges, and then cross out 
all resulted isolated vertices, the resulting graph will be connected. 
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Proof. Suppose the assertion to be false. Denote by V ne (G) the set of all 
non-end vertices of G. Then for every vertex s € V ne (G) the result of crossing it 
out together with all its edges, followed by crossing out all isolated vertices, is a 
graph G s with several components: G s = G Sj i U G s ,2 U • ■ ■ U G S: m s with m a > 2. 
Let e m i n (s) be the minimum number of edges in G s .i, i = 1,2, ...,m„. Define 
so £ V ne (G) to be a minimizer of the function e m ; n (■) . Fix a component G° S(j 
with e m i n (so) edges. Note that G — G° is connected. Take any vertex v\ 6 G2 
which is a neighbor of s . Clearly, si G Kie (G). Let us cross it out from G. By 
assumption, G Sl also has several components, the one containing G — G° o is not 
the smallest one, hence the smallest one lies in G° o and has therefore less edges 
than e m i n (sq), which brings us to contradiction. ■ 

Proof of the Theorem We argue by induction on the number of edges, 
and we use the notation of the proof of the Lemma ^ above. The case of one 
edge is clear. By Lemma ^| we can find a non-end vertex s 6 V ne (G), such 
that the cross-out graph G s remains connected; since a (G) < a (G s ) + 1 and 
e (G s ) < e (G) — 2, we get the first four statements. 

To prove (|2^) we first prove the upper bound. We have: < \ + 2eTG) • 

Note that always e (G) > v (G) — 1, therefore < \ + 2 ( v (G)-i) > wnne 

e (G) > u (G) and ^1^1 < i + wnenever there is a cycle in G. To prove 

that the left hand side of ( p6f ) is less than or equal to the right hand side, the 
only case left is when G is a tree with v (G) odd. This is done by induction on 
v. For v = 3 the statement is clear. If v > 3, fix a vertex so £ V (G) and take 
an end vertex si most distant from so ■ The only vertex S2 adjoint to it has 
the graph G S2 connected. We have a (G) < a (G S2 ) + 1, e (G) = v (G) — 1 and 
e (G S2 ) = v (G S2 ) — 1 (since both G and G S2 are trees). If w (G S2 ) = v (G) — 2, 
then by the induction hypothesis 

a (G) < a (G S2 ) + 1 < — - — + + 1 



2 2«(G S2 ) 



2 V «(G S2 ) 



2 V «(G)y 2 \ v(G) 

If v (G S2 ) < v (G) — 3, then similarly 

a(G)<a(G„) + l<— + 
= i(«(G. B ) + 2) 



2 \ v (G) J 2 V « (G) 
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Since the estimate (g6|) is achieved on the line graph (for v even) and on the 
circle graph (for v odd), the proof is complete. ■ 

Acknowledgments. We would like to thank V. Kleptsyn, whose proof 
of Theorem |s| is reproduced above, V. Dolotin for his help with Fig. 1 and 
L.Bassalygo for valuable remarks. 
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